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Abstract 

For two-loop two-point diagrams with arbitrary masses, an algorithm 
to derive the asymptotic expansion at large external momentum squared 
is constructed. By using a general theorem on asymptotic expansions of 
Feynman diagrams, the coefficients of the expansion are calculated ana- 
lytically. For some two-loop diagrams occurring in the Standard Model, 
comparison with results of numerical integration shows that our expansion 
works well in the region above the highest physical threshold. 
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1 Introduction 



The high precision that has been achieved in high energy e + e~ experiments, espe- 
cially at LEP, makes them sensitive to radiative corrections in the Standard Model 
at the one and two loop level. However, some complicated theoretical problems are 
connected with the evaluation of two-loop corrections involving massive internal par- 
ticles (heavy quarks, W and Z bosons, Higgs particles) whose masses cannot be 
neglected in the region of interest. In particular, much attention has been paid to 
the evaluation of massive two-loop self-energy diagrams. In some well known special 
cases, such as the QED correction to the photon self-energy, they can be evaluated 
exactly Jl|, § , and the result can be expressed in terms of trilogarithms. The problem 
becomes essentially more difficult when all the internal particles of the diagram are 
massive. Exact results for the corresponding Feynman integrals are not known, and 
we need to look for other approaches to calculate the contributions of such diagrams. 

One of the possible ways is by numerical integration. For this purpose, it is 
convenient to use a two-fold integral representation || obtained for the so-called 
"master" two-loop two-point integral. Another way is by constructing asymptotic 
expansions of such diagrams in different regions, the coefficients of the expansions 
being calculated analytically. For example, in ref. @] we constructed an expansion 
for the general case of the two-loop self-energy diagram that works when the external 
momentum squared (k 2 ) is below the first physical threshold. We showed that, unless 
we are very close to the threshold, only a few terms of the expansion are needed to 
obtain highly accurate results. 

In the present paper we shall consider another situation, when k 2 is larger than the 
highest threshold of the diagram. This case is more complicated, because the corre- 
sponding expansion is not a usual Taylor expansion, but also contains logarithms and 
squared logarithms of {—k 2 ) (in four dimensions) yielding an imaginary part when 
the momentum is time-like. The procedure of calculating this diagram in the form of 
a series in inverse powers of the external momentum (plus logarithms) resembles the 
standard procedure of analytic continuation of the hypergeometric function (some ex- 
plicit examples of such a procedure connected with Feynman integrals were presented 
e.g. in ref. ||). To obtain this expansion, we shall apply a general mathematical the- 
orem on asymptotic expansions of Feynman integrals in the limit of large momenta. 
This theorem holds at least in the case when the external momenta are not restricted 
to a mass shell. However, in our case of a two-point massive diagram, it is valid for 
any values of the external momentum. All the coefficients of the expansion have a 
natural interpretation in terms of Feynman integrals and are analytically calculable 
in the case we consider. Expansions of this kind were presented in refs. || [7|, |^. || 
and rigorously proved in ref. |10| (see also ref. |Tl| for a review). 

The remainder of this paper is organized as follows. In section 2 we explain how 
the general theorem can be applied to construct the large momentum expansion for 
the general mass case of two-loop self-energies. We present an algorithm by which 
any term of this expansion can be calculated analytically. In section 3 we discuss the 
analytical results obtained for the coefficients. Section 4 is devoted to a comparison of 
the asymptotic expansion with the results of a numerical integration using the method 
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of ref . . Section 5 contains our conclusions. Finally, a number of formulae we need 
for massless two-point integrals and for massive vacuum integrals are collected in 
appendices A and B, respectively. 



2 Constructing the asymptotic expansion 

We shall consider the "master" two-loop self-energy diagram presented in Fig. la. The 
corresponding Feynman integral can be defined as 



/ / d n p d n q 

J{{vi}]{m}]k)= I I D ^ D v2 D v, D v, D ^ (i) 



where Di = (pf — mf + iO) are massive denominators (for brevity, we shall omit the 
"causal" zO's below), Ui are the powers of these denominators, n = 4 — 2e is the 
space-time dimension (in the framework of dimensional regularization fT3, 13]), and 



the momenta of the lines pi are constructed from the external momentum k and the 
loop integration momenta p and q (with due account of the momentum conservation 
in all vertices). Two-loop diagrams with three or four internal lines correspond to 
the cases when some of the z/s in eq. (|l]) are equal to zero. Moreover, in the case 
of integer z/s the two- loop diagram with another topology (see Fig. lb) can also be 
reduced to integrals fll]), if we use the decomposition of first and fourth denominators 
(see e.g. in 0). So, if we deal with integer powers of denominators, it is in general 
sufficient to consider the integrals (|l]) with different powers z/j. 

Now let us introduce the result that the asymptotic expansion theorem gives for 
our case (the general case of this theorem can be found in ref. []TT|): 



00 J2 J ?h ° 2{mi}i{«}«Jy> ( 2 ) 

k 2 — >oo ^ 

where T is the main graph (see Fig. la), 7 are subgraphs involved in the asymptotic 
expansion (see below), r/7 is the reduced graph obtained from T by shrinking the 
subgraph 7 to a single point, J 7 denotes the dimensionally-regularized Feynman in- 
tegral corresponding to a graph 7 (for example, Jr is given by ([!])), Ti mi yr qi ^ is the 
operator of Taylor expansion of the integrand in masses and "small" momenta qi 




(a) (b) 
Figure 1: Two-loop self-energy diagrams 
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(that are "external" for the given subgraph 7, but do not contain the "large" ex- 
ternal momentum k), and the symbol "o" means that the resulting polynomial in 
these momenta should be inserted into the numerator of the integrand of Jr/ 7 . It is 
implied that the operator T acts on the integrands before the loop integrations are 
performed. 

In our case, the sum (0) goes over all subgraphs 7 that become one-particle ir- 
reducible when we connect the two vertices with external momentum k by a line. 
In other words, we should consider all possibilities to "distribute" the external mo- 
mentum k among and for each arrangement the subgraph 7 will coincide with 
the subset of lines involving k, while lines without k should be removed. For our 
graph T (see Fig. la), all possible subgraphs 7 (there are five different types of them) 
are presented in Fig. 2 (dotted lines correspond to the lines that do not belong to 
7). 




(7 = r); 



Figure 2: The subgraphs 7 contributing to the large k 2 expansion 

The reduced graphs r/7 correspond to the dotted lines and can be obtained by 
shrinking all solid lines to a point. In such a way, we obtain that for the second and 
third type (see Fig. 2) J r / 7 corresponds to a massive tadpole, for the fourth type we 
obtain a product of two massive tadpoles, while for the fifth type we get a two-loop 
massive vacuum integral (with three internal lines). 

The Taylor expansion operator T expands the denominators of the integrand in 
the following way: 



T 



OO 

\ " 



[p 2 



m 



2\3 



(3) 



1 



m;q 



[(k-q) 2 -m 2 } u [k 2 ] 1 



where 



j=0 J' 

I> + j) 

r(iz) 



u)j ( 2(kq) — q 2 + m 



2\3 



k 2 



(4) 
(5) 
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is the Pochhammer symbol . In fact, on the r.h.s. of (|j) it is understood that the real 
expansion goes over the total power of m and q (so, for example, the terms q 2 , m 2 and 
(kq) 2 /k 2 should be considered together). The reason is that when we evaluate the 
massive vacuum integrals Jr/ 7 with these momenta q in the numerator, the powers 
of q will be transformed into powers of the masses. If we have several denominators 
to expand, then we should also collect all terms with the same total power of masses 
and "small" momenta gj. 

Now we are able to consider what integrals correspond to different terms of the 
asymptotic expansion (see Fig. 2). 

1. In this case, 7 = Y. All denominators of (HD should be expanded in masses, 

T {nH} J{{v l } ] {m l } ] k)= £ ■ ■ ■ fok* (mi)* ■ ■ ■ (m 5 )* J({u i + jj; {0}; k). 

ji,...,j 5 =0 31 J5- 

(6) 

Note that if we consider the case V\ — . . . — v§ — 1, the first term of the expansion 
(H) (with ji — . . . — j 5 — 0) gives the well-known result: — 6C(3)-7r 4 //c 2 . Two-loop 
massless integrals with higher integer powers of denominators occurring on the r.h.s. 
of (|6|) can be evaluated by use of the integration-by-parts technique [14j] (see Appendix 
A). 

2. Let us consider only the first contribution of the second type (see Fig. 2), when 
7 is obtained from T by removing line 1. Then we get 

d> f d n q 



[ P 2 -m 2 } Ul m2 -~y [ g 2- m f[(j,-g)2- m 2pp-p)2- m |]^p- g )2_ m 2p' 

(7) 

After expanding the integrand of the g-integral in masses and p, we obtain products of 
massless one-loop integrals (see Appendix A) and massive tadpoles with numerators 
that can be calculated by use of eq. (|34D (see Appendix B). Other contributions of 
the second type can be evaluated in the same way. 

3. In the case when the central line is removed from 7, we obtain: 

d n p 



[p 2 — m 2 }" 3 

X Tm 



mi ' m2 ' m4 ' m5;p i [(k+p~q) 2 -m 2 } Ul [(k-q) 2 -m 2 } U2 [(p-q) 2 -m 2 } Ui [q 2 -m 2 } U5 ' V J 

After expansion, we obtain integrals of the same type as in the previous case. 

4. There are no loop integrations in the subgraph 7, and we get for the first 
contribution of the fourth type: 

d> d n q 



\p 2 — m 2 ] Ul \q 2 —m 2v '' 



1 



XTm2,m3 ' m4M li(k-q) 2 -m 2 p[(k-p-q) 2 -m 2 pi(k-p) 2 -m 2 p) ' (9) 

As a result, we obtain products of two one- loop tadpoles with numerators (also for 
the second contribution), which can be evaluated by use of eq. (|35|) (see Appendix 
B). 
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5. In this case we obtain a non-trivial two-loop vacuum integral. For example, 
the first contribution of the fifth type gives: 



d n p d n q 



lp 2 -m 2 } vl [q 2 -m 2 } U2 [(p-q) 2 ~™ 2 ^ 3 



rrin 



xTm4 '"™ \{{k-pY-mT[{k-q) 2 -mT ) ' (10) 

Expanding the denominators, we obtain two-loop vacuum integrals with numerators, 
that can be evaluated by use of eq. fl36[) presented in Appendix B (the same for the 
second contribution of the fifth type). Note, that here we obtain the same two-loop 
vacuum integrals as those involved in the small- k 2 expansion (see |4[]). In particular, 
for unit powers of denominators the dependence on masses can be expressed in terms 
of dilogarithms (see (|2|)-(|33|)). 

So, the total asymptotic expansion is the sum of all terms presented, and we see 
that all integrals corresponding to the coefficients of the expansion can be evaluated. 



3 Analytical results 

In principle, eqs. (|^)-([10|) presented in the previous section (together with the for- 
mulae of Appendices A and B) enable one to construct analytical expressions for 
the coefficients of the large-fc 2 expansion. However, in the general case of unequal 
masses the higher-order coefficients become rather cumbersome. To calculate these 
coefficients, we used the REDUCE system [TjJ. The algorithm constructed is applica- 
ble to integrals with arbitrary values of masses, space-time dimension and (integer) 
powers of denominators. If we are interested in the result near n = 4, we perform an 
expansion in e = (4 — n)/2 to get the divergent and finite parts of the coefficients. 

One of the most important examples is the "master" two-loop diagram (presented 
in Fig. la) in the case v\ = . . . = v§ = 1. In this case, the result should be finite 
as n — > 4 (and it is a non-trivial check of the algorithm that all the divergent 
contributions from separate terms of fl6|)— QTU) cancel in this sum !). A rigorous proof 



of the finiteness of the expansion was given in || (it was based on the so-called 
iT-operation [[Uj). 
Let us define 

J(l,...,l;mi,...,m 5 ;fc)| n=4 = J(m u . . . ,m 5 ;fc) 

= JL M {m x ,...,mxk) = -^^^, (11) 

where the coefficient functions M.j include the powers of masses and the logarithms of 
masses and momentum squared. It is easy to see that the only integral contributing 
to Mq is J(°)(l,l,l,l,l) in dD (see eq. (|26D in Appendix A). So we get the obvious 
result that the expansion starts from 

= 6C(3). (12) 
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The Aii term already includes contributions of all terms (|6|)-(|T0|) (with the ex- 
ception of ® that begins to contribute starting from M.?)] this yields 

Ml = ^/ 111 2f_4u 41n f_4 N )_ ln 4 111 !4 +6 i 

2 I V mf I \ mf) mf mf j 



+ | analogous terms with 
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( 2 In* + 4 In (-^) - In 4 In 4 - In 4 In 4 } 

2 1 \ m| ) V m\ J m| m\ m\ m\ J 

+- {F[m\, m 2 2 , ml) + F{m\, m\, ml)) , (13) 

where the symmetric function F is defined by fl30l)-(|33|) (note that F has the dimen- 
sion of mass squared). We see that M.\ contains also ln(— k 2 ) and In (—k 2 ) terms, 
which means that our expansion (at n = 4) is not a usual Taylor expansion. In fact, 
the highest power of ln(— k 2 ) is connected with the highest order of pole in e (double 
pole) that can occur in the two-loop integrals involved. Note that for positive k 2 the 
sign of the imaginary part produced by these logarithms is defined by the "causal" 
iO-prescription, 

ln(-A; 2 - zO) = ln(fc 2 ) - m ik 2 > 0). (14) 

Let us also present the result for the next term of the expansion, 

M 2 = 4W(-4U 4 l„(-4)-21n4l„4 + 7} 
8 [ \ mf) \ mf J mf mf J 

■i 1 



+ | analogous terms with m\, m\, m\ 



m\{ , 9 / k 2 \ ( k 2 \ , m? , mi , m\ , ml 
+ { -2 In 2 ? - 2 In ? + In — i In — § + In — | In — | + 6 

4 [ y m| J \ m| / m| m| m| m| 

/ 2 2 i 2 2\ 

- -Am x m 2 + m 4 m 5 ) 

+ + h > f-^ + 4!n (-4) + 4,n (- 

\ mf) \ mi) 



•(■ 


fc 2 \ 
m 2 J 


+ In 2 (- 


fc 2 \ 
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mf 


mf 


m| 


m% 



- In -f In -§ - In -f In -f + i 

mf mf m| m| 

+ | analogous term with m 2 m 2 | 



+ I analogous term with m^m 2 } 

m 2 ml [ n , 9 / A; 2 \ , / k 2 \ , m? m 2 , , m5 , m= 1 
+ -^-^2 In 2 f -2 In ? - In —| In — 5: - In — | In — | - 8 f 

2 1 \ m 3 / \ m l / m 3 m 3 m 3 m 3 J 

+ I analogous terms with m^m 2 , m\m\, m 2 777 2 j 

+ — I (^m 2 + m\ — m 2 + 2m 2 + 2m^j F(m 2 , m 2 ,, m|) 

+ (2m 2 + 2m 2 , - m 2 + m 2 + m 2 ) F(m 2 , m 2 , m 2 )} . (15) 

Higher contributions are more cumbersome, and we do not present them here (but 
we are going to use them below, when comparing our expansion with the results of 
numerical integration). By use of the REDUCE system ||15|| , for the general massive 
case of the integral ([TT[) (when all five masses are arbitrary) we obtained analytical 
results for the coefficient functions up to A4q. 

There are also some other possibilities to check the correctness of our results (in 
addition to cancellation of 1/e poles). For example, in ref. analogous results were 
presented for some special cases when some of the masses are zero while others are 
equal (see also ref. |L7j where these results were generalised to the case of arbitrary 
space-time dimension n). In these cases, we compared our results for the coefficients 
(at n — 4) and found complete agreement. 



4 Numerical results 

In this section, we will continue to focus on the "master" two-loop integral correspond- 
ing to Fig. la. In general, it has two two-particle thresholds, at k 2 = (mi + 7714) 2 and 
k 2 = (m 2 + rris) 2 , and two three-particle thresholds at k 2 = (mi + m 3 + m 5 ) 2 and 
k 2 = (m 2 + m 3 + m.4) 2 . The asymptotic expansion theorem quoted in (Ef) provides a 



series of approximations to this integral of the form (see eq. (11)): 



N KA ■ 



% (k 2 Y 



such that the remainder behaves like 

M - M {N) = O {(k 2 )- N - 1 ln 2 (-A; 2 )) (17) 

as k 2 — > 00. 

Strictly speaking, this does not imply convergence of the series for any fixed value 
of k 2 , but from experience with special cases where exact analytical results are known, 
one would expect the series to converge when \k 2 \ is larger than the highest threshold. 
In order to see whether this is true in the general case, and whether the asymptotic 
expansion can be used as a practical means of calculating these integrals, we made 
some numerical comparisons for two physical examples. 
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The first example is: 



J(M w ,M w ,M z ,m b ,m b ;k), (18) 

where the M^y, M z and m b denote the masses of the P^-boson, the Z-boson and the 
bottom quark. This integral is a contribution to the top quark self-energy. In this 
case, both the two-particle thresholds coincide. The three-particle thresholds also 
coincide with each other. 

The other example occurs in the photon, the Higgs and the Z-boson self-energies: 

J(m t ,mt,Mz,mt,m t ;k). (19) 

Here m t is the top quark mass. As in the first example, there are only two distinct 
thresholds. We calculated ([18]) and (|i~9f ) numerically using the method of ref . . The 
values we took for the masses were: 

M z = 91 GeV, M w = 80 GeV, m t = 140 GeV, m b = 5 GeV. (20) 

The results are displayed in Figs. ^ and || In each plot, the first threshold is exactly on 
the left edge and the position of the highest threshold is indicated by a dashed vertical 
line. The dotted horizontal line shows the lowest order asymptotic approximation 
A^ -* = AAq = 6^(3) = 7.21234 . . . (the imaginary part is zero to this order). The 
curves labeled 1, 2, . . . show the approximations A^ 1 -*, AA^ 2 \ . . . defined by eq. (|I~6|). 
The results of the numerical integration program are shown as crosses. 

At large values of k 2 , say ten times the highest threshold or higher, Ai^ already 
agrees with the numerical results to within 0.01%, which is the order of magnitude 
of the error in the numerical results. 

The most interesting region is immediately above the highest threshold, where we 
would still expect the expansion to converge, but more slowly than at large k 2 . In fact, 
in our first example (|18j), A^ 4 -* is still accurate to within 1% on the threshold itself. 



In the second example fll9[) , the convergence near the threshold is less rapid. This 
region is shown in more detail in Fig. |5], where the left edge of the plots corresponds 
to the highest threshold. On the threshold itself, the error in Al (6) is about 6% m 
the real part), but it drops down to less than 1% by the time k 2 is a factor of 1.2 
above the threshold. 



Note that in |THJ radiative corrections to the on-shell H — > bb squared amplitude 



(represented by a diagram like Fig. la) were examined, and it was pointed out that 
the mass correction A^i essentially improved the description above the last threshold. 

Finally, when we go below the highest threshold, we see large deviations between 
the asymptotic approximations and the numerical values, and the series ceases to 
converge. From this we see that the asymptotic expansion can only be applied to the 
region above the highest threshold, as expected. 

5 Conclusions 

Thus, in the present paper we considered an algorithm to construct an asymptotic 
expansion of two-loop self-energy diagrams (see Fig.l) for large values of k 2 . This 
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algorithm can be applied to the general case of different masses, and integer powers of 
denominators. We used the asymptotic expansion theorem fl2|) to derive the different 
contributions to the expansion (see (|6|)-(|T0|)). In any order, the coefficients of the 
expansion can be expressed in terms of known one- and two-loop propagator-type 
massless integrals, and one- and two-loop massive vacuum integrals. 

By use of the REDUCE system, we wrote a program that automatically generates 
analytical expressions for the coefficients of the asymptotic expansion. Then we 
considered the "master" two-loop diagram with different masses (Fig. la), and we 
obtained expressions for the terms of the asymptotic expansion of the integral ([TTJ) 
up to l/(k 2 ) 7 terms (with logarithms). The first three coefficient functions A^o,i,2 
are given by equations (|12D, ([13|) , ([15]). For some concrete diagrams occurring in 
the Standard Model, we made a numerical comparison with the results of a two-fold 
numerical integration based on the algorithm of ||. We found that in the region 
above the highest physical threshold of the diagram, our expansion converges well, 
and in the region not very close to the threshold it is sufficient to take only a few 
terms of the expansion to provide a good accuracy. 

Thus, the present paper (together with ||J]) solves the problem of constructing 
asymptotic expansions of massive two-loop self-energy diagrams in the limits of small 
and large k 2 values (when we are either below the lowest physical threshold or above 
the highest one). The analytical description of the behaviour between the thresholds 
still remains a problem. Note, however, that the general asymptotic expansion theo- 
rem can also be used for other limits in order to calculate the considered diagram in 
situations when some of the masses are large and the others are small. 

A procedure analogous to the one considered here can be also applied to three- 
point two-loop vertex diagrams, when all external lines are above the corresponding 
thresholds (and also in some other cases). Note that some two- loop massless integrals 
with off-shell external momenta (needed for such an expansion) were calculated in 



ref. [ 19 1 . For numerical comparison, a parametric integral representation obtained in 



ref. p0| can be used. 



Acknowledgements 

We are grateful to F.A. Berends and W.L. van Neerven for their assistance and helpful 
discussions. We are also grateful to D.J. Broadhurst, R. Buchert, A.L. Kataev and 
D. Kreimer for useful discussions, and to H. Kuijf and S.C. van der Marck for writing 
the program that calculates the integral (|TT|) numerically. One of the authors (A.D.) is 
grateful to Instituut-Lorentz for their hospitality and to the Nederlandse organisatie 
voor Wetenschappelijk Onderzoek (NWO) for financial support. Another author 
(V.S.) is thankful to the Institute for Theoretical Physics of Gottingen University for 
kind hospitality and to the Heraeus Foundation for support. 



10 



Appendix A: Massless integrals 

Here we present some formulae needed for evaluation of one- and two-loop massless 
integrals occurring in the paper. These results are well known, and we write them 
only for completeness. 

The massless one-loop integral is 



J(0)( 



•l-n„n/ 2/ ,„2W2-^ 2 r (?V 2 ~ jj£Wg ~ ^("l + "2 ~ Tl/2) 



where n = 4 — 2e is the space-time dimension. 

Massless two- loop integrals are defined by (see ([!])) 

J (0) (i/i, u 2 , u 3 , 1/4, fs) = J{{vi}] {0}; fc) 



(22) 



By use of symmetry properties of the diagram in Fig. la, and the integration by parts 
formula 



J {0) (vi, V2, V3, "4, vs) = (fi + 2z/ 3 + - n) 1 

X fa J (0) (z/l + 1, V 2 - 1, ^3, ^4, ^5) - J^\vi + 1, V 2 , ^3 ~ 1, ^4, t> 5 ) 



,^2,^3,^4 + 1,^5 - 1) - J(°)(l/ a ,V 2) I/3 -l,I/ 4 + 1,^5) ] } , (23) 



these integrals (with positive integer z/'s) can be reduced to the following "boundary" 

J (0) K^,0,^ 5 ) =I®{yi,v 4 ) / {0) K^); (24) 

J^fa, v 2 , v 3 , i/ 4 , 0) = I®(v 2 , v 3 ) J(«0(i^ + i/ 2 + 1/3 - n/2, z/ 4 ). (25) 

For example, 

j(°)(l l l l l,l,l) = ~6C(3) + 0(e) (26) 

(note that £(3) does not occur in the divergent and finite parts of any other integrals 
with positive z/'s). 



Appendix B: Massive vacuum integrals 

The result for one-loop massive tadpole integral is well known Jl2| : 

Two-loop vacuum massive integrals 

f f d n p d n q 

I(vi,V2,V3;m 1 ,m 2 ,m 3 ) = / / y- 2l — j-= 2 -— r ^ gr— (28) 

J J [p 2 — mf\ Ul [q 2 — m 2 \ V2 [{p — q) 2 — m 3 \ U3 
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have been studied e.g. in ref. 0] (some special cases were also considered in ref. |2~l| ). 
For example, for unit powers of denominators we get 



1(1, 1, l;m 1 ,m 2 ,m 3 ) 



71 



4-2e-p2 



T 2 (l+e) (l + 3£ + 7£ 2 



X 



j — -^-^(ml + ml + m\) H — (m^ lnm 2 + m\ mm 2 , + m\ lnm 2 ) 



1 r 
2 



2 1 2 2, 2 1 2 2, 2 1 2 2 

m l in m l + m 2 in m 2 + m 3 in m 3 



+ (mf + — mj) In ml In m 2 . + (ml — m 2 2 + m 2 ) In ml In mg 
+ (— ml + m 2 , + mj) lnm 2 , lnm 2 + F(ml, ml, m 2 ) j + 0(e) 



(29) 



where the function F is symmetric with respect to m\, mf>, tttJ, and is defined by (see 
I) 

F (ml, ml, ml) = m\\ 2 (x,y)§ (x,y) (30) 



with 



x 



m% 



ms 



ni 



2 ■ 



A 2 (x, y) = (l - x - t/) 2 - 4xy, 



(31) 
(32) 



$(x,7/) = - 21n( ^ In ^ j-lnxlny 



-2Li 2 



1 + x — y — y 



If the largest mass is mi or m2, we should choose this mass as the dimensionless- 
making parameter in (30) and (^l|) (instead of 777.3). m the region where A 2 < 0, the 
function (RDT) can be represented in terms of Clausen's functions (see M). To obtain 



results for integrals fl28|) with higher integer powers of denominators, a recursive pro- 
cedure based on integration by parts [l4j can be applied (see in f4|). It is interesting 
that the same function ( |33| ) also occurs when we evaluate massless triangle diagrams 
(see e.g. in |22|; in this case x and y are dimensionless ratios of external momenta 
squared) . 

We also need to evaluate massive integrals with numerators. To do this, in the 
one-loop case we used the following explicit formula: 



d n p 



[p 2 



m 



2]u 



[2(Mf 1 [2(Mf 2 



N!+n 2 - even 



iVi! iV 2 ! 



(n/2) (M 



(N!+N 2 )/2 



E 

2jl+j S =Nl 
I 2j 2 +j 3 =V 2 



(kiy\kiy 2 [2(h k 2 )} 

ii! h s - J3! 



[P 2 



777 



2]u 



(P 



2\(N 1 +N 2 )/2 



(34) 



12 



and the integral on the l.h.s. vanishes if (Ni + iV 2 ) is odd. The sum in braces goes 
over all non-negative integers ji,j2,33 obeying two conditions: 2j 1 + j 3 = N% and 
2j2 + 33 = ^2 (so, in fact it is a one- fold finite sum). We also use a standard notation 
for the Pochhammer symbol fl5|). 

When we have a two-loop integral with a numerator and two denominators (cor- 
responding to a product of two one- loop tadpoles), the following formula can be 
derived: 



d n p d n q 



[2(kp)] N i [2{kq)]»> [2{pq)\ 



[P 2 



m\} vi [q 2 



mi 



\U2 



N 2 



]N 3 



N!+N 3 - even 
N2+N3- even 



iVi! N 2 \ (k 



2 x(Vi+V 2 )/2 



{n/2) {Ni 



+N 2 )/2 



£ 

2jl+33=Ni 
I 2j 2 +h=N 2 



(N3+33) 



Ji!j2!j 3 ! ((iV 3 +j 3 )/2)!(n/2) (iV3+i3)/2 
d n p d n q 



[p 2 — m 2 } Ul [q 2 — m 2 ] V2 



{p 2f 1+ N 3 )/2 {q2) (N 2+ N 3)/ ^ (35) 



and the integral on the l.h.s. of ( p5[) is equal to zero if (JVi + N 3 ) or (iV 2 + N 3 ) is odd. 
The sum in braces is of the same structure as in (J34f) . 

We also need an analogous formula for the two-loop integral with three denomi- 
nators: 



d n p d n q 



[p 2 - m 2 ] Ul [q 2 - m%\ V2 \{p - q) 2 - m^ 3 



[2(kp)}^ [2(kq)} 



N 2 



^+N 2 - even 



2 x(jVi+AT 2 )/2 



N2) (Ni+N2)/2 2ji+js=m tititi 

2j 2 +j 3 =N 2 

x / / 21 r 9 ^fu v> 21- (pTilT [2(p?)] i8 , (36) 

J J [p 2 -m 2 ]^[q 2 -ml]^[(p - q) 2 -mj] u * v ' K ' L v n v 7 

and this integral is also equal to zero if (N\ + N 2 ) is odd. 

In all formulae (|34|), (|35D and (^), the remaining momenta in the numerators 
on the r.h.s. can be expressed in terms of denominators, and we arrive at the result 
expressed in terms of integrals (1271) and without numerators. 
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Figure captions 



Figure 3: The real and imaginary parts of J(M W , M w , M z , m b , m&; k) for 
k 2 > (M w + m b ) 2 . The second threshold is at (M w + M z + m h ) 2 . 



Figure 4: The real and imaginary parts of J(m t , m t , M z , m t , m t ; k) for k 2 > Am 2 . The 
second threshold is at (2m t + M z ) 2 - 



Figure 5: The real and imaginary parts of J(m t ,m t , Mz,m t ,m t ;k) for 
k 2 > (2m t + M z ) 2 . 
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